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Abstract. 

We study the possibility of nanoscale phase separation in manganites in the 
framework of the double exchange model. The homogeneous canted state of this model 
is proved to be unstable toward the formation of small ferromagnetic droplets inside an 
antiferromagnetic insulating matrix. For the ferromagnetic polaronic state we analyze 
the quantum effects related to the tails of electronic wave function and a possibility of 
electron hopping in the antiferromagnetic background. We find that these effects lead 
to the formation of the threshold for the polaronic state. 



1. Introduction 

Manganites, the Mn-based magnetic oxide materials such as LaMn03, have been known 
for more than 50 years. Jonker and van Santen pQ, and in more details Wollan and 
Koehler investigated a rich magnetic structure of Ca-doped La 1 _ ;E Ca a ;Mn03 and 
conductivity in these materials. Specifically, a strong correlation between magnetic and 
transport properties in manganites was observed. These materials exhibit a metal-like 
resistance in a ferromagnetic phase and an insulating behavior in an antiferromagnetic 
phase. 

While the variety of magnetic structures in manganites was explained by 
Goodenough j3] based on the theory of semicovalent exchange, the correlation between 
transport and magnetic properties was for a first time qualitatively explained by Zener 
He suggested that the conduction electrons travel in Lai_ x Ca x Mn03 through 
the Mn 4+ ions and each ion carries a fixed magnetic moment which is strongly 
ferromagnetically coupled with the spin of a carrier by a generalized Hund rule. Since 
a spin of a conduction electron should be aligned parallel to the local spin, in the 
classical picture a conduction electron cannot move in antiferromagnetically ordered 
environment. Anderson and Hasegawa 5 j solved a problem of two local spins and 
one conduction electron, providing a strong mathematical background to Zener ideas. 
They found that if the bare hopping amplitude was t, than a hopping amplitude of 
an electron moving between the two local spins is t e ff = tcos9, where 9 is an angle 
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between the directions of the local spins. Thus, the kinetic energy of conduction electrons 
directly depends on an angle between the sublattice moments. Moreover, the conduction 
electrons tend to align ferromagnetically the local spins, which surround them. Hence, 
a competition between ferromagnetic coupling via a conduction electron (so called 
double exchange mechanism) and an antiferromagnetic coupling via superexchange of 
two neighboring local spins takes place. De Gennes [S| suggested that this competition 
results in a homogeneous canted state, namely an angle 9 is uniform and monotonically 
changes from ir (collinear antiferromagnetic order) to (collinear ferromagnetic order) 
with increasing carriers concentration x. Based on these considerations de Gennes 
plotted the first phase diagram of the manganites. 

Soon after this work, Nagaev [7J |H] improved de Gennes results by considering 
quantum fluctuations associated with the local spins. Nagaev proved that electron 
can move even in antiferromagnetically ordered phase with a small hopping amplitude 
tj \J2S + 1. Nagaev, Kasuya, and Mott also proposed [HI HOI UJ] that for small 
electron concentrations it is more favorable for conduction electrons to form a self- 
trapped state (ferromagnetic polaron or ferron) in antiferromagnetic matrix by creating 
a small ferromagnetic bubble, rather than to form a homogeneous canted state in the 
whole sample. Thus, it was one of the first hints for nanoscale phase separation in a 
double-exchange model made by that time. 

Recent growth of interest in manganites was initiated in 1993 by the discovery of 
the colossal magnetoresistance (CMR) effect in doped LaMnOs. The CMR phenomena 
implies a drastic decrease of resistivity in manganites in the presence of a magnetic field 
|12[ IT3*]. Soon after discovery of CMR in manganites a phase diagram of La 1 _ 2 ,Ca ;c Mn03 
was revised |14j . It was found that a variety of phases appears in manganites in addition 
to those predicted by de Gennes. 

There are plenty of experimental and theoretical studies of manganites nowadays. 
They were initiated first of all by the potential technological applications of colossal 
magnetoresistance phenomena and also by the interesting physics of strong correlations, 
which manifests itself in these materials. In particular, the interaction of spin, charge, 
and orbital degrees of freedom in manganites as well as their rich phase diagrams draw 
much attention of theorists and experimentalists in recent years. 

The important question that has to be answered is about a leading mechanism 
responsible for CMR in the optimum doping region (x ~ 0.3). Some authors argue 
that CMR could be explained in framework of the double-exchange mechanism alone 
|15j . others claim that it is necessary to take into account a lattice interaction 
(Jahn- Teller polarons), some insist that percolation-type arguments could explain CMR 
|17j . However, as it was pointed out by Dagotto et al jTH] and Arovas and Guinea 
[T9~] both analytical and numerical calculations in various models related to manganites 
exhibit a strong tendency toward phase separation in the wide range of temperatures 
and concentrations. Thus, it is believed that CMR phenomena could be understood as 
a competition and co-existence of different phases in manganites as well as a phenomena 
related to the proximity of the optimum doping region to various phase-separation 
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thresholds. Note that at higher concentrations close to half-filling there appears another 
threshold of phase separation in the system corresponding again to the formation of 
ferromagnetic droplets, but now in a charged ordered insulating matrix [_U] . 

One of the authors j2Ij demonstrated that the double-exchange model at low doping 
is unstable toward phase separation and the energy of a homogeneous canted state is 
higher than an energy of a self-trapped state corresponding to a ferromagnetic polaron. 
Hence a legitimate question arises whether the stability of a polaronic state is preserved 
when quantum fluctuations of spins and tails of wave function of conduction electrons 
are taken into account. 

2. Basic theoretical model 

The general chemical formula of the most popular class of manganites is Lnx_2.Aj.MnO3, 
where Ln is a trivalent cation Ln 3+ (La, Pr, Nb, Sm, . . .), and A 2+ is a divalent cation 
(Ca, Sr, Ba, Mg, . . .). Oxygen is in a 2 ~ state, and the relative fraction of Mn 4+ and 
Mn 3+ is controlled by a chemical doping x. This class of manganites has the perovskite 
structure. In the cubic lattice environment, the five-fold degenerate 3<i-orbitals of Mn- 
ions are split into three lower-energy levels (d xy , d yz , and d zx ), usually referred to as t2 g , 
and two higher-energy e g states (d x 2_ y 2 and d 3z 2_ r 2). The t<i g levels with three electrons 
form a state with a local spin S = 3/2, whereas delocalized e g states contain an electron 
or are unoccupied depending on the chemical doping x. The e g states are further split by 
the static Jahn - Teller effect and for the simplicity we will treat here only the lowest e g 
state, assuming a Jahn- Teller gap to be large enough and neglecting any orbital effects 
in our consideration. 

The simplest theoretical model suggested for the explanation of the properties of 
manganites is the ferromagnetic Kondo lattice model (s — d model): 



The first term in equation represents a strong on-site Hund's ferromagnetic coupling 
(Jh > 0) between the local spin 5 = 3/2 and the spin of a conduction electron. In 
real manganites, the Hund's interaction Jh is of the order of 1 eV. The second term in 
equation Q is the kinetic energy of the conduction electrons. The projection operator 
P corresponds to the case of singly occupied e g orbitals (a strong Hund's interaction 
prevents two conduction electrons with different spin projections from occupying the 
same site). Note that a strong electron-lattice interaction significantly reduces the 
effective width W of the conduction band (W = 2zt) resulting in a rather small hopping 
amplitude t « 0.3 eV. The third term in equation (0) is a weak antiferromagnetic 
coupling between local spins on neighboring sites, with Jff ~ 0.001 eV. In equation ((TJ, 
symbols mean the summation over z nearest neighbor sites. 

In the case of a strong on-site Hund's coupling (Jh _> W J//) the model 
described by the first two terms of the Hamiltonian (P) is referred to as the double 
exchange model. Note that if all local spins are ferromagnetically aligned, the 




(1) 
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conduction electrons will move freely in their surrounding. Thus, model (JTJ describes the 
competition between the direct antiferromagnetic coupling of local spins and the double 
exchange via conduction electrons, which tends to order local spins ferromagnetically. 
In the strong-coupling limit, Hamiltonian Q can be simplified: 

h = - J2 *(%K +a ; + J ff s2 E cos (%) ( 2 ) 



(i,J) 



where af and a,j are creation and annihilation operators of spinless fermions (conduction 
electrons whose spins a are aligned parallel to the local spins), £(%) is an effective 
hopping amplitude, and 9y is an angle between sublattice moments, as we already 
discussed. The hopping amplitude in the case of classical spins (S ^> 1) reads: 

t(9) = tcos(0/2) (3) 

In Nagaev's quantum approach the local spins at empty sites have the maximum 
projection, +S, on the magnetization vector of the corresponding sublattice. At 
occupied sites, however, local spin S and spin cr of a conduction electron form a state 
with the total spin S + 1/2, but with two possible projections of it S ±1/2. So there are 
two effective bands in quantum canting corresponding to the two different projections 
of the total spin. Their bandwidths read [HI I2T]: 

t 



t±(0) 



^2S + 1 + S 2 cos 2 (#/2) ± S cos(fl/2) 



(4) 



25 + 1 

The quantum hopping amplitude drastically differs from the classical de Gennes one. In 
contrast to the de Gennes picture, an electron can still move through antiferromagnetic 
matrix creating a state with S* ot — S + 1/2 at one site and a state with S^ ot = S — 1/2 
at a neighboring site (a string-like motion introduced by Zaanen and Oles 22 J: 



lot 



s + 



s. 



tut 



s 



s. 



tot 



s 



Hence, there are two equal hopping amplitudes in the case of electron motion through 
antiferromagnetic background: t + = t_ = t/y2S + 1. On the other hand, for 
ferromagnetic ordering one gets from equation t+ — t and £_ = t/(2S + 1). The 
proportionality of an effective bandwidth in the case of quantum canting to l/VS is 
just an implication of its quantum nature. 



3. Homogeneous canted state 

An energy of the classical de Gennes state with an account taken for the 
antiferromagnetic interaction between the local spins reads: 

E = -ztx cos(9/2) + ~zJf f S 2 cos(9), (5) 

where z is the number of nearest neighbors, and x is a carriers concentration. The first 
term in this equation is the gain in the kinetic energy, and the second term is the loss 
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X c l = 

3 



(7) 



in the energy of antiferromagnetic interaction between local spins. Minimization of the 
energy (jSj) with respect to the parameter cos(0/2) yields: 

Thus, we have a transition from a collinear antiferromagnetic state for x = to a 
collinear ferromagnetic state for = 2JffS 2 /t. For < x < x C 4, the canting angle 

(9 ^ 7r), and a homogeneous canted state takes place. 

Previously, various homogeneous states were considered taking into account the 
quantum hopping amplitudes (El El]. In contrast to the classical collinear 
antiferromagnetic state remains energetically favorable up to the critical value of the 
carrier concentration x c \, which is given by: 

UffS 2 i I 3 
zt + 1_ 

Thus, in a quantum case, the canted state occurs for x > x c \, whereas in the classical 
case the canted state appears for arbitrarily low doping levels. At higher doping levels 
(x > x c i), the two-band quantum canted state arises, namely conduction electrons are 
in the two bands with total spin projections \S^ ot = S + 1/2) and \S^ ot = S — 1/2). 
However at x > x C 2 ~ (27/2)x c \ the bottom of the second band lies above the chemical 
potential level and a one-band state of quantum canting becomes favorable [21]. Finally, 
at much higher carriers concentration (x > x C 3 = AJffS 2 /t\/2S + 1), a transition to the 
classical canted state of de Gennes (0) occurs. Note that for x > x C 4 = 2JffS 2 /t the 
canted state transforms into a collinear ferromagnetic state with the angle 8 = 0. 

To test the stability of the homogeneous state the electronic compressibility was 
calculated [2] according to the standard formula = d 2 E/dx 2 . At very low 
concentrations x < x c \ the compressibility is positive and a homogeneous collinear 
antiferromagnetic state corresponds to at least a local minimum of energy. However in 
the range of concentrations x c \ < x < x c ^ the compressibly k < and the canted state 
is unstable. For example, we can calculate the compressibility of the classical de Gennes 
canted state and get: 

, zt 2 

Hence compressibility of this state is negative, i.e. de Gennes classical canted state 
is also unstable. The negative sign of the compressibility indicates the instability 
of homogeneous state toward the phase separation. The simplest case of nanoscale 
phase separation corresponds to the formation of small ferromagnetic polarons inside 
an antiferromagnetic matrix. This state was considered in Ref. |211 I2H] and it was 
shown that a polaronic state is more favorable energetically than all the homogeneous 
states in the total range of concentrations < x < x^. Note that magnetic polarons, in 
this case, correspond to the electron in the self-trapped ferromagnetic state of a finite 
radius inside an antiferromagnetic insulating matrix. 
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4. Polaronic state 



As we already discussed, in the case of the classical hopping amplitude (JBJ) a conduction 
electron may be self-trapped and form ferromagnetic droplets (magnetic polarons) inside 
an antiferromagnetic matrix. The simplest assumption is to consider that the boundary 
between the ferromagnetic region and the antiferromagnetic matrix is abrupt without an 
extended region of inhomogeneous canting. Then an energy of a polaronic state reads: 



E 



-tx z 



R 2 



1 T Q2 4vr (R 



2 zJ ff S 



4tt (R 
1 - -z-x - 
3 \ a 



(8) 



In equation (JSJ), R is a radius of a polaron and a is a lattice constant. The first 
term in equation (JSJ) describes the kinetic energy gain due to the formation of a 
ferromagnetic region. The corrections to this term proportional to ta 2 /R 2 correspond to 
the localization energy of a conduction electron inside a ferromagnetic droplet of radius 
R. The second term in equation (JSJ) is a loss in the Heisenberg antiferromagnetic energy 
of local spins inside the droplet. Finally, the third term describes the energy of an 
antiferromagnetic interaction between local spins in a region outside the ferromagnetic 
polarons. The polaron radius is obtained from the condition of energy minimization 
dE/alR = 0. So we have the following expressions for an energy and a polaron radius: 



E 



pol 



— ZtX H 7T tX 



Tit 



2/5 



n zJ ff S 



R 



pol 



Tit 



1/5 



(9) 



(10) 



2zJffS\ 

Note that in this case the transition from a polaronic to a ferromagnetic state 
occurs when ferromagnetic polarons start to overlap. The critical concentration for 
ferromagnetic transition reads: 

3 /fl x3_ 3 ( 2zJffS^ 3/5 m) 

In the case 



Xc5 



47T \RJ 47T V Kt J 

Now, let us consider the quantum corrections to these results, 
corresponding to the quantum hopping amplitude (J2J), a conduction electron can move 
through antiferromagnetic background with a heavy mass m* ~ yS, so it is interesting 
to study for this case a problem concerning the stability of a magnetic polaron. Since 
the analysis of the discrete model (J2J) is rather complicated, we consider the continuum 
limit assuming that a radius of a polaron is much larger than the lattice constant a 
(further on, we put a = 1). A total energy (J2J) can now be written in the following form 



E = - J [z\i)\ 2 + i)*/\i)]t{6)dV + zJf f S 2 J cos 2 (9)alV, 
t{9) 



t 



2S+1 



^2S + l + S 2 cos 2 (#/2) + S cos(fl/2) 



(12) 
(13) 




Figure 1. The electron wave function ip{r) and the canting angle 9{r) in the quantum 
case (solid lines) and the classical case (dashed lines) . 



As one can see from equation (|12p. the total energy will lie between the two 
limiting values corresponding respectively to the motion of the conduction electron via 
ferromagnetic or antiferromagnetic backgrounds: 

zt 



Efm = —zt < E < E AFM 



V2S + 1 

Since the electron wave function should be normalized J \ip\ 2 dV = 1, we minimize the 
functional F = E — t/3 J \ip\ 2 dV with respect to parameters 9 and if), where j3 is a 
Lagrange multiplier. The corresponding Euler-Lagrange equations have the following 
form: 



[2zip + AiP] t{6) + A [t(6)ip] - 2f3ty = (14) 

sin(0/2) = (15) 



[z\^ + rM) d ^ m - 2zJ ff S 2 cos(e/2) 

We solve these two coupled differential equations by the following iterative procedure 
[2*1] : (a) we choose a trial function for the canting angle 9(r); (b) we solve the first 
differential equation ()14j) to obtain an electron wave function ip(r); (c) using an obtained 
value for ip(r) we solve equation (|15|) to get a canting angle function 9(r); (d) we return 
to the step (a) until our iteration process converges. 

Functions ip(r) and 9(r) obtained by the numerical solution of equations ()14j) and 
(|15J) are shown in figure Q for a broad range of the values of parameter a = t/JffS 2 . 
Note that exact numerical solution which takes into account both the effects of quantum 
canting and tales of the wave function coincides with the classical Nagaev-Mott solution 
for a — ► oo. One can see that the magnetic polaron represents a very good localized 
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Figure 2. The ground state energy corresponding to polaron formation in the 
quantum case (open circles) and the classical case (solid circles). For local spin S = 3/2 
ferromagnetic and antiferromagnetic states correspond, respectively, to Efm = — 6t 
and Eafm = —2>t. 



object, and the transition region from the ferromagnetic ordering (6 = 0) to an 
antiferromagnetic matrix (9 = tt) is narrow enough. Nevertheless, a polaronic state can 
disappear at relatively small value of the parameter a c ~ 75. Indeed, as one can see from 
figure |2l there is a transition from the polaronic state to a collinear antiferromagnetic 
state at small values of parameter a < a c . For this case, the total energy of a magnetic 
polaronic state is equal to the bottom of the band for electron motion through the 
antiferromagnetic background, and as a result for a < a c an electron can move freely 
through the antiferromagnetic matrix. Note that to get a more precise value of a c we 
should solve a variational problem for the functional F on the discrete lattice since for 
small values of a a continuous approximation is not accurate enough. The work along 
these lines is in progress now. 



5. Conclusion 



The tendency toward phase-separation is very strong in the double-exchange model. We 
have shown that in the wide range of parameter a = t/JffS 2 & conduction electron forms 
a self-trapped state. In this state, an electron is localized in the ferromagnetic droplet of 
finite radius embedded in the antiferromagnetic matrix. This construction seems rather 
natural in the de Gennes classical approximation of the double exchange, where hopping 
amplitude t e ff = tcos(6/2) and electron cannot move through an antiferromagnetic 
background since t e ff = for 9 = tt. However, we have proved that even in the quantum 
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case, when a conduction electron can travel slowly through antiferromagnetic matrix 
(since t e ff = t/y/2S + 1 for 9 = n), the polaronic state remains well defined and stable. 
Our approach to the one electron problem in an antiferromagnetic matrix corresponds 
to a very small doping region in real materials, where the concentration of charge 
carriers (e.g. holes in La x Cai_ x Mn03) is low. However we believe that even for higher 
concentrations the polaronic picture remains qualitatively correct. Moreover, at very 
low concentrations magnetic polarons should be localized on impurity sites [213 whereas 
at concentrations higher than the Mott threshold they are depinned from the impurities. 
Preliminary estimates show that the Mott threshold in manganites corresponds to 
xm ~ 0.01 — 0.04, which is significantly lower than the critical concentration x c § ~ 0.15 
for the overlap of magnetic polarons. 

Our model of polarons embedded in ferromagnetic matrix allows one also to 
calculate a magnetoresistance and a noise spectrum of manganites, if we suggest that an 
electron transport takes place via the hopping of electrons from one polaron to another, 
while a polaron itself is immobile. These calculations were carried out in Ref. [2E1 Elj 
and found experimental support in the recent paper of Babushkina et al [21] . 

If we proceed now to the experimental confirmation of small-scale phase separated 
picture we should mention that there already exist many evidences in favor of nanoscale 
phase separation in low and moderatly doped manganites. The confirmation of 
an inhomogeneous state in manganites comes from the nuclear magnetic resonance 
experiments of Allodi et al \2Q^iQ, where the two different hyperfine lines corresponding 
to ferromagnetic and antiferromagnetic regions were observed. Experiments on neutron 
scattering of Biotteau et al [31 J support the idea of small ferromagnetic droplets 
embedded in an antiferromagnetic or a canted matrix. And if we turn to transport 
properties, a very natural picture of electron percolation, which is in agreement with 
our model, was experimentally confirmed by Babushkina group in Ref. [32J. 

Thus the experimental and theoretical picture strongly confirms a phase-separated 
state in manganites in the region of low doping. Moreover, we believe that this picture 
remains qualitatively correct for the concentrations optimal for the CMR effect in the 
high temperature region T > Tq (Tc is a Curie temperature), where the ferromagnetic 
fluctuations of the short range (the temperature polarons) are present j2HllSS|- Hence, a 
combination of very intuitive picture of polarons and the ideas of the percolation theory 
could provide a correct description for the behavior of manganites in the wide range of 
temperatures and carrier concentrations. 
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